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Abstract 



We have studied the properties of A = 54 and A = 64 isobars at tempera- 
tures T < 2 MeV via Monte Carlo shell model calculations with two different 
residual interactions. In accord with empirical indications, we find that the 
symmetry energy coefficient, b sjm , is independent of temperature to within 
0.6 MeV for T < 1 MeV. This is in contrast to a recent suggestion of a 2.5 MeV 
increase of 6 sym for this temperature, which would have significantly altered 
the supernova explosion scenario. 
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There appears to be a general consensus on the basic scenario for presupernova collapse. 
An iron core is produced in the center of a massive star in the final stage of hydrostatic 
nucleosynthesis. When the iron core exceeds the Chandrasekhar mass limit, it can no longer 
be supported by electron degeneracy pressure and so begins to collapse. Subsequent pho- 
todisintegration of nuclei and electron capture on nuclei (later on free protons) reduce the 
pressure and energy of the core and accelerate its collapse, in this way raising the central 
density and temperature. After neutrino trapping, the inner part of the core (mass densi- 
ties p > 10 12 g/cm 3 ) collapses as a homologous unit whose size is given by the appropriate 
Chandrasekhar mass limit, Mch ~ 5.76(Y" e trap ) 2 M , where Y^ rap is the electron-to-nucleon 
ratio of the trapped material. The collapse halts when the homologous core exceeds nuclear 
matter density (p ^ 10 14 g/cm 3 ). A shock wave, formed at the inner core's surface, then 
travels outward and eventually explodes the star. 

Despite the appeal of this scenario, the actual mechanism of a supernova explosion is still 
controversial. In the direct mechanism the shock wave is strong enough not only to stop the 
collapse, but also to explode the outer shells. In the delayed mechanism the shock energy 
is first stored in the core mantle (the layer of infalling material in the iron core outside the 
homologous region), but it is recovered shortly thereafter by the outstreaming neutrinos. A 
detailed description of the physics of a supernova can be found in Ref. |TJ. 

The size of the core mantle is a major determinant of the supernova mechanism. It 
has been argued recently that the mantle size might be significantly smaller than generally 
calculated in supernova models because of an overlooked temperature dependence of the 
electron capture process @. The detailed reasoning is as follows ||: When the nucleus is 
described by a static mean field model, dynamical effects associated with collective surface 
vibrations are conventionally embodied in an effective nucleon mass, m* (as are spatial non- 
localities in the mean field). Donati et al. || studied the coupling of the mean field single- 
particle levels to the collective surface vibrations within the quasiparticle random phase 
approximation and found that m* decreases with increasing temperature for T < 2 MeV, 
the temperatures relevant for the presupernova collapse. They attributed this behavior to a 
reduction of collectivity at low excitation energies. In the Fermi gas model, this decrease of 
m* induces an increase in the symmetry energy contribution to the nuclear binding energy 

E sjm {T) = 6 sym (T)- — , (1) 

where N, Z, and A are the neutron, proton, and mass numbers of the nucleus, respectively. 
Quantitatively, Donati et al. estimate that b sym (T) increases by about 2.5 MeV as T increases 
from to 1 MeV (6 sym (0) « 28 MeV |§). 

Importantly, a larger value of E gym (T) would hinder electron capture in the presupernova 
environment, and thus increase Y" g tra P. This, in turn, would reduce the size of the core mantle 
so that the shock wave would need significantly less energy to stop the collapse and explode 
the star. In the supernova scenario outlined by Donati et al, the temperature dependence 
of 6 sym increases the shock energy by about 0.5 foe, roughly half of the observed explosion 
energy of SN87a [§. 

A more definitive test of the proposed temperature dependence of the nuclear symmetry 
energy (and hence of the presupernova electron capture process) necessarily involves a proper 
account of those nuclear degrees of freedom beyond the mean field (e.g., pairing correlations) 
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that are relevant for nuclei under presupernova conditions (e.g., /p-shell nuclei at T < 
2 MeV). The model of choice for this task is the interacting shell model, which provides a 
complete nuclear spectrum and so allows a statistical description using a canonical ensemble; 
both quantal and thermal fluctuations (within the model space chosen) are fully included. 

The recent development of the Monte Carlo shell model || has made feasible complete 
(Ohu)) calculations of /p-shell nuclei at finite temperature; a pilot study of the thermal 
properties of the nucleus 54 Fe has been presented in Ref. [El. To explore the temperature 
dependence of the symmetry energy, we have used this method to study the thermal proper- 
ties of two isobar chains (A = 54 and 64) in the mass region important for the presupernova 
collapse; this includes two of the three nuclei studied in Ref. || ( 64 Zn and 64 Ni). Technical 
considerations restrict our Monte Carlo shell model calculations to the N = Z and even-even 
members of the chains. The details of our calculation parallel those of the 54 Fe study in 
Ref. H and the Monte Carlo shell model is described in detail in Refs. |5|,|7]]. For the A = 54 
isobar chain (o 4 Co, o 4 Fe, ^ 4 Cr, and | 4 Ti) we have adopted the Brown- Richter interaction 0. 



For the A 



and ° 4 Fe), 



for which the Brown-Richter force is 
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64 isobars («Ge, «Zn, ^Ni, ^ 38 , 
expected to be unreliable (it was constructed from properties of nuclei at the beginning 
of the fp-shell), we have used the KB3 residual interaction | I~0|] . Further, to explore the 
sensitivity of our results to the two-body interaction, we have also studied the A = 54 chain 
with the KB3 interaction. All of our calculations omit the Coulomb interaction and assume 
complete isospin symmetry. 
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Figure 1. Calculated internal energies, U, for the A = 54 and A = 64 isobar chains at various 
temperatures using the KB3 interaction. Results for A = 54 with the Brown-Richter interaction 
are virtually identical. The energies are given relative to the doubly-magic nucleus 40 Ca. The 
experimental energies have been corrected for the Coulomb repulsion as described in the text and, 
in the case of A = 64, have been increased by 11 MeV to account for the systematic underbinding 
of the Hamiltonian for these isobars. 



In Fig. 1 we display our results for U, the internal energy (expectation value of the 
Hamiltonian in the canonical ensemble) for the A = 54 isobars and A = 64 isobars as a 
function of (N — Z) 2 at temperatures relevant for the presupernova. For even-even nuclei 
our calculation at T = 0.5 MeV has almost cooled to the ground state, as the first excited 
states in these nuclei have excitation energies E x ^ 1 MeV. For the N = Z nucleus 54 Co, the 
T = 0.5 MeV calculation represents a mixture of the ground state and the first excited state 
at 200 keV. If we correct the experimental binding energies by the semi-empirical Coulomb 
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energy, E c = 0.717Z 2 /-4 1/3 (1 - 1.69/A 2 / 3 ) MeV Q, we find that our calculation with either 



interaction accurately reproduces the A = 54 binding energies. For the A = 64 nuclei, the 
low-temperature calculations are systematically overbound by some 11 MeV, although the 
variation along the chain is well reproduced. A better description of the nuclei 64 Ni and 64 Fe 
(with 16 and 18 valence neutrons in the /p-shell) would plausibly require inclusion of the 
gg/2 orbital. We observe that, at a given temperature, U is linearly proportional to (N — Z) 2 
for the even-even nuclei of an isobar chain. This behavior is, of course, expected from semi- 
empirical parametrizations of the binding energy (e.g., the Bethe-Weizsacker formula) and 
almost exclusively reflects the symmetry energy contribution to the binding energies (1). 

For A = 54, our calculation yields 6 sym (0) = 21 ±1.0 MeV and 22±0.75 MeV for the KB3 
and Brown- Richter interactions, respectively. These values are in very good agreement with 
experiment, 6 sym = 21 MeV. For A = 64, we find b sym = 23 ± 0.7 MeV, slightly larger than 
the experimental value of 21.2 MeV. We note that all of these values are significantly smaller 
than the symmetry energy coefficient adopted in presupernova studies (6 sym = 30-32 MeV 

00)- 



Our calculation does not confirm the temperature dependence of the symmetry energy 
proposed in Ref. ||. As is obvious from Fig. 1, the slope of the linear variation of U with 
(N—Z) 2 for the even-even nuclei is independent of temperature, indicating that 6 sym does not 
change. To be more quantitative, an increase of 6 sym by 2.5 MeV between T = and 1 MeV, 
as proposed in Ref. ||, corresponds to an increase in the energy splitting between 54 Fe and 
54 Ti of A = 4.4 MeV (using the experimental value for 6 sym ), while our calculation yields 
A = 0.4±0.9 MeV, implying that b sym changes by less than 0.7 MeV. In fact, our calculation 
shows that all properties of the even-even nuclei studied here are essentially constant for 
T < 1 MeV. These include the quadrupole strength, whose variation with temperature has 
been cited as the source for the proposed temperature-dependence of the symmetry energy 
in Ref. §. 

The low-temperature constancy of the properties of even-even nuclei is expected in view 
of their sparse spectrum of low-lying excited states and relatively high excitation energies of 
their first excited states. To confirm this, we have calculated the nuclear thermodynamics 
from the low-lying experimental levels of the nuclei in the two isobar chains; enough of these 
levels are known to give reliable results up to T = 0.7 MeV. The partition function is given 
by 

Z(T)=J2(2J+l)e- Ej/T , (2) 
J 

where the free energy is given by F(T) = — T In Z(T), and the internal energy is U = 
—T 2 d(F/T)/dT. We then calculated the difference in internal energies between isobars a 
and b as a function of the temperature, and expressed that in terms of 6 sym through Eq. (1). 

For the A = 54 chain, we used all known levels in 54 Fe and 54 Cr up to 10 MeV of 
excitation energy; levels with unknown spin were assigned a J-value randomly chosen from 
the standard empirical spin distribution function |jl3"l. Results for U are shown in the left- 



hand panel of Fig. 2, expressed in terms of A6 sym , the change in b sym relative to its value at 
T = 0. There is structure in the curve due to the peculiarities of the 54 Fe and 54 Cr spectra, 
although it is smaller than the precision of our Monte Carlo calculations. We find that \b sym \ 
changes by 0.3 MeV at most as T varies from to 0.7 MeV. A similar calculation in the 
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A = 64 chain using 64 Zn and 64 Ni gives a change in |6 sy m| of 0.6 MeV as T varies from to 
0.6 MeV. In contrast, Donati et al. claim that b sym increases by 1.22 (1.53) MeV for 64 Ni 
( 64 Zn) as T varies from to 0.7 MeV.Q Note that although we have defined the symmetry 
energy coefficient in terms of U, rather than F, the right-hand panel of Fig. 2 shows that 
the symmetry coefficient of F never varies by more than 0.5 MeV for T < 0.7 MeV. It is 
also interesting to note that, depending on the low-lying spectra, 6 sym may either increase 
or decrease as T increases. 
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Figure 2. Temperature dependence of the symmetry energy coefficients for the internal energy, U, 
and free energy, F, calculated from the empirical spectra as described in the text. 54 Fe and 54 Cr 
were compared to calculate the A = 54 case, while 64 Zn and 64 Ni were used for A = 64. The limited 
knowledge of levels at high excitation energy makes these cases unreliable for T > 0.7 MeV. 



In contrast to even-even nuclei, odd-odd nuclei (with iV not equal to Z) have a rather 
rich spectrum at low excitation energies. Thus, U is expected to increase more strongly for 
these nuclei at temperatures T < 1 MeV. From a simple Fermi gas picture one expects an 
increase of aT 2 »s 7 MeV, with the empirical level density parameter a = A/8 MeV -1 [13 



The increase of U with temperature for odd-odd nuclei changes the kinematics of electron 
capture, making it easier for odd-odd nuclei to capture electrons. 

As is demonstrated in detail in Ref. || for the case of 54 Fe, our Monte Carlo shell 
model calculations predict the disappearance of pairing correlations between like nucleons 
at T « 1.3 MeV, while np-correlations persist to higher temperatures. We observe the 
same behavior for the nuclei studied here. At T « 1.3 MeV we find a significant change in 
the thermal properties of the nuclei. For example, the moment of inertia and the (orbital) 
Ml strength increase drastically. This change is also apparent from the increase of U with 
temperature, as shown in Fig. 1. Monitored by the number of BCS-like pairs in the nuclei, 



1 We have calculated the change of the symmetry energy coefficient predicted by Donati et al. 
using the exponential fit formula given in Table I of Ref. Q. We note, however, that this formula 
significantly underestimates the change in 6 S y m at T — 1 MeV compared to the QRPA results given 
in the same table. 
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we observe that the change is related to the disappearance of the pp- and wz-correlations, 
which, in turn, dominate the nuclear properties at low temperatures. The vanishing of the 
pairing energy in even-even nuclei at T > 1.3 MeV has a significant effect on the kinematics 
of the electron capture on nuclei in the late stage of presupernova collapse. 

Another novel feature is apparent in our calculations, although it has no direct relation 
to the presupernova situation. It is well known that odd-odd N = Z nuclei are slightly less 
bound relative to the even-even nuclei in an isobar chain because of the lack of like-particle 
pairing of the last proton and neutron. However, there is apparently a strong correlation 
between these unpaired nucleons in N = Z nuclei that energetically favors them over the 
other odd-odd nuclei in the chain. For both interactions, we observe in the A = 54 isobar 
chain that 54 Co becomes energetically favored relative to the even-even nuclei, after the 
disappearance of the like-pair correlations in the latter. This observation is in agreement 
with the finding in |J that isoscalar np- correlations persist to higher temperatures than the 
like-pair correlations. Apparently these correlations still contribute to the internal energy 
of 54 Co at T > 1.3 MeV. 

In conclusion, Monte Carlo shell model calculations show clearly that pairing correlations 
play a decisive role in nuclear properties at low temperatures. For all (N = Z and even- 
even) /p-shell nuclei studied, we find that the nuclear properties are essentially constant 
for T < 1 MeV, as expected from the spectrum of these nuclei. In particular, we find that 
the symmetry energy coefficient is constant to within some 0.6 MeV for T < 1 MeV. Thus, 
our calculation does not confirm speculations based on a mean-field approximation that the 
symmetry energy increases by 2.5 MeV in the temperature interval T = 0-1 MeV. Our 
calculations illustrate that nucleon-nucleon correlations beyond the mean-field level must be 
taken into account. A detailed account of this work will be given elsewhere. 

This work was supported in part by the National Science Foundation, Grants No. 
PHY90-13248 and PHY91-15574. We thank the Concurrent Supercomputing Consortium 
for a grant of DELTA and PARAGON time. 
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